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ABSTRACT 



From 1979-1996 the Green Bank Interferometer was used by the Naval Research Laboratory 
to monitor the flux density from 146 compact radio sources at frequencies near 2 and 8 GHz. 
We filter the "light curves" to separate intrinsic variations on times of a year or more from more 
rapid interstellar scintilation (ISS) on times of 5-50 d. Whereas the intrinsic variation at 2 GHz 
is similar to that at 8 GHz (though diminished in amplitude) , the ISS variation is much stronger 
at 2 than at 8 GHz. We characterize the ISS variation by an rms amplitude and a timescale 
and examine the statistics of these parameters for the 121 sources with significant ISS at 2 GHz. 
We model the scintillations using the NE2001 Galactic electron model assuming the sources are 
brightness- limited. 

We find the observed rms amplitude to be in general agreement with the model, provided 
that the compact components of the sources have about 50% of their flux density in a component 
with maximum brightness temperatures 10 11 -10 12 K. Thus our results are consistent with cm- 
wavelength VLBI studies of compact AGNs, in that the maximum brightness temperatures found 
are consistent with the inverse synchrotron limit at 3 x 10 11 K, boosted in jet configurations 
by Doppler factors up to about 20. The average of the observed 2 GHz ISS timescales is in 
reasonable agreement with the model at Galactic latitudes above about 10°. At lower latitudes 
the observed timescales are too fast, suggesting that the transverse plasma velocity increases 
more than expected beyond about 1 kpc. 

Subject headings: ISM: general — scattering — plasmas — galaxies: active — galaxies: quasars — 
galaxies: BL Lacertae objects 
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1. INTRODUCTION 

Variations in the radio flux density of extra- 
galactic sources have been observed for many 
decades; see early observations by Dent (1965). 
Via the classic light travel time argument, vari- 
ability timcscales have been used to infer the 
physical size of the radio emitting regions, and 
the associated apparent brightness temperature. 
Kellerman & Pauliny-Toth (1968) reviewed the 
early work on centimeter wavelength variations in 
terms of radiation from (non-relativistic) expand- 
ing synchrotron sources. In 1981 the same authors 
reviewed the topic in the light of very long base- 
line interferometry (VLBI) observations of appar- 
ent superluminal expansion. Since that time a 
standard picture has emerged of relativistic jets 
from active galactic nuclei (AGNs), which emit 
Doppler-boosted radio beams. If the beam points 
close to the Earth, these are seen as very compact 
sources with fiat radio spectra. The relativistic 
analysis of such beams successfully explains the 
sporadic brightcnings which emerge from a cen- 
tral core and travel out "supcrluminally." In some 
models the brightcnings are thought to be shocks 
and in others physical structures moving with the 
jet. 

This now canonical jet model has typical 
Doppler factors in the range 1-20. Flux variations 
are converted to a "variability brightness temper- 
ature," calculated by the light travel time argu- 
ment; this increases as the cube of the Doppler fac- 
tor over the ~ 10 12 K limit due to inverse Comp- 
ton scattering (e.g., Blandford & Konigl 1979) 
or over the <~ 3 x 10 11 K equi-partition limit as 
proposed by Readhead (1994). Typical year-long 
flux variations at centimeter wavelengths give ap- 
parent brightness 10 14 -10 15 K, which can be ac- 
comodated within the expected range of Doppler 
factors. However, variations of flux at frequencies 
below about 1 GHz proved harder to fit in to the 
jet model. For a source of a given flux density and 
angular size, the associated brightness tempera- 
ture scales as the inverse square of the frequency 
Thus the observation of substantial variations over 
a year at a frequency of 0.4 GHz (e.g., Fanti et 
al. 1981) pushed the apparent brightness to over 
10 16 K and stretched the limits of the relativistic 
beaming models, by requiring higher Doppler fac- 
tors to explain the low-frequency variations (LFV) 



than to explain the superluminal expansions. 

While the discovery of pulsars confirmed that 
radio sources can vary intrinsically on very short 
timescales, their flux variations over minutes to 
hours were quickly found to be due to scintil- 
lation caused by inhomogenities in the interstel- 
lar medium (ISM). As reviewed by Rickett (1977) 
these phenomena were found to be consistent with 
a very wide range of length scales in the ion- 
ized ISM, which causes strong scintillations (rms 
> mean). Dcnnison & Condon (1981) looked for 
and set upper limits to such scintillations from a 
variety of compact extragalactic sources including 
several LFV sources. The absence of interstellar 
scintillation (ISS) must be due to quenching by 
the angular extent of the emitting regions, from 
which they inferred lower limits to the angular 
size, which were larger than those inferred from 
an intrinsic explanation of LFV. 

This paradox was resolved by the recognition of 
a further regime of ISS. At meter wavelengths pul- 
sars show flux variations over days to months as 
well as over minutes. Sieber (1982) demonstrated 
that the slower variations were also a propaga- 
tion phenomenon, which Rickett, Coles & Bour- 
gois (1984) identified as predicted by strong scin- 
tillation theory. This is now referred to as re- 
fractive interstellar scintillation (RISS) to distin- 
guish it from the diffractive scintillation (DISS) 
over minutes to hours. Since the angular diame- 
ter that quenches RISS is substantially larger than 
for DISS, these authors proposed that RISS could 
also account for LFV. The suggestion that LFV 
was due to irregular interstellar propagation had 
earlier been made by Shapirovskaya (1978), who 
proposed that refraction in discrete large scale ion- 
ized structures was responsible. Several authors 
(e.g., Rickett 1986; Spangler et al. 1993) have 
since demonstrated that RISS does indeed pro- 
vide a satisfactory explanation for LFV, based on 
a power-law model for the spectrum in the inter- 
stellar plasma density. 

A further variability phenomenon at centimeter 
wavelengths was discovered by Heeschen (1984). 
In careful daily measurements he found a "flick- 
ering" at a level of a few percent over 1-2 days 
in the radio flux from many flat spectrum sources. 
The sources which showed the greatest amplitudes 
were followed up in extensive monitoring cam- 
paigns, particularly, by the Bonn group. Quir- 
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renbach et al. (1992) presented the results of daily 
monitoring of a complete sample of flat spectrum 
sources at 6 and 11 cm. They found a low level 
of what they called intraday variability (IDV) in 
most of the sources observed, which also had com- 
pact VLBI structure. They concluded that this 
low level of variabilty seen in almost all such 
sources was RISS. However, in addition they ar- 
gued that about 25% of flat spectrum sources also 
showed intrinsic variations that can reach am- 
plitudes as high as 20%. An intrinsic variation 
of 0.2 Jy on a timescale of 12 hr implies a vari- 
ability brightness over 10 19 K, implying a Dopplcr 
factor over 200. Thus for those sources said to be 
showing intrinsic IDV, flux variations were inter- 
preted in terms extreme brightness emission that 
strains the limits of the relativistic jet models. 

In fact there should be a continuum of ISS be- 
havior between LFV and IDV. The wavelength 
dependence is determined by a combination of the 
steep increase with wavelength for the character- 
istic scale of RISS and the wavelength scaling of 
effective source diameter. For most extragalactic 
sources their angular size partially quenches the 
scintillations and the effective scintillation ampli- 
tude and time scale are determined by the compact 
source size, which scales approximately as wave- 
length (for brightness temperature limited emis- 
sion). In an analysis of the flux variations of 
quasar 1741—038 over a wide range of frequen- 
cies, Hjcllming & Narayan (1986) identified ISS at 
a frequency of 1.49 GHz and slower intrinsic vari- 
ations at 15 and 22 GHz. Dcnnison ct al. (1987) 
showed that RISS was the dominant cause of the 
2.7 GHz variation of two sources at low Galactic 
latitude, which were in the Green Bank Interfer- 
ometer monitoring program. 

A detailed analysis of one of the large amplitude 
IDV sources (0917+624) by Rickett et al. (1995) 
demonstrated that ISS could successfully account 
for the IDV with a source brightness of about 
6 x 10 12 K, providing that the scale height of the 
scattering medium was reduced to 200 pc, from the 
500 pc scale height of the Taylor & Cordes (1993) 
model for the Galactic distribution of electrons. 
This was one of the sources for which 20% intrin- 
sic IDV had been argued earlier. 

An epoch of extremely rapid variability was re- 
ported for quasar 0405—385 by Kedziora-Chudczer 
et al. (1997). They observed large amplitude vari- 



ations on timescales as short as one hour at 4.8 
and 8.4 GHz. They considered the implied vari- 
ability brightness 10 21 K to be unreasonable and 
argued for an ISS explanation. Rickett et al. 
(2002) reported detailed analysis of the IDV in 
the polarization of this source, which they found 
to be consistent with a maximum brightness tem- 
perature of 2 x 10 13 K. 

Two other quasars have also shown substantial 
variation over times as short as an hour (intra-hour 
variation - IHV), namely J1819+385 (Dennett- 
Thorpe & de Bruyn, 2000, 2002 & 2003) and 
B1257-326 (Bignall ct al. 2003). The measure- 
ments of a time delay in the light curves recorded 
on trans-continental baselines and the observa- 
tion of annual modulations in their characteris- 
tic timescales give irrefutable evidence for scin- 
tillation as the cause of the variations for these 
two quasars. It is also important to note that the 
three IHV sources mentioned here require that the 
distance to the scattering medium be only tens 
of parsecs rather than the 500-2000 pc distances 
involved in most pulsar ISS. A consensus has fi- 
nally emerged that the IDV (and IHV) phenom- 
ena are due to ISS. The strongest argument in 
support of intrinsic IDV comes for one particular 
quasar (0716+714), for which 6 cm IDV was ob- 
served to have some correlation in behavior with 
optical IDV observed simultaneously (Wagner et 
al. 1996). However, the argument is weakened 
by the fact that the optical and radio variations 
were not correlated in detail and there was signif- 
icant decorrelation across the radio wavelengths 
observed (see observations by Quirrenbach et al., 
2000). 

Variations over months and years have been 
well established at shorter wavelengths, as for ex- 
ample reported by Stevens et al. (1994) who mon- 
itored a group of blazars for several years at 22 to 
375 GHz. They found reasonable agreement with 
shock models for the many bursts they observed. 
They also observed faster variations at lower lev- 
els, but the brightness temperature deduced by 
assuming intrinsic variations are not extreme at 
their high frequencies. 

Most variability studies are based on relatively 
small samples of sources, often observed at differ- 
ent telescopes using different observing method- 
ologies. Lazio et al. (2001, henceforth L01) pre- 
sented the light curves and structure functions 
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for a sample of 146 compact radio sources moni- 
tored at the Green Bank interferometer (GBI) over 
long durations (3-15 yr) with frequent flux density 
measurements and a common observing methodol- 
ogy. They concluded that the rms variation in flux 
density at both frequencies was dominated by slow 
large amplitude changes, due mostly to intrinsic 
variability on timescales from a few months to a 
few years. They also found more rapid variations 
on timescales of 5-50 days, which they character- 
ized as being ISS. 

This paper reports further analysis of the GBI 
flux monitoring program. We shall show that the 
variations at 2.25 GHz do indeed form a contin- 
uum with the ISS phenomena in LFV and IDV. 
The plan of this paper is the following. In §2 we 
summarize the GBI monitoring program. In §3 we 
describe our analysis. In §4 we describe a model 
for the ISS, which we compare with the observa- 
tions in §5. in §6 we discuss the results for a few 
sources of particular interest and we give our con- 
clusions in §7. 

2. Observations 

Flux density observations with the GBI were 
made nearly daily from 1979-1987 on 46 radio 
sources, as described by Fiedler et al. (1987b, 
henceforth F87) and Waltman et al. (1991, hence- 
forth W91). The sample of sources was expanded 
to 146 and the frequency of observations reduced 
to once every 2 days during 1988-1996 (L01). Here 
we summarize the relevant details of the Navy- 
NRAO GBI monitoring program; for complete de- 
tails see F87, W91, and L01. 

Sources were selected for monitoring on the ba- 
sis of their known variability in flux density and/or 
compact structure. Flux densities were obtained 
at two frequencies in S-band and X-band. Un- 
til 1989 August the frequencies were 2.7 GHz (S 
band) and 8.1 GHz (X band) after which they were 
changed to 2.25 GHz (S band) and 8.3 GHz (X 
band) when new cryogenic receivers were installed. 

The observations were made on a baseline 
of 2.4 km and the correlated flux density esti- 
mated from 15 min. integrations of 35 MHz band- 
width summing right and left circular polariza- 
tions. Each source was normally observed near 
the same hour angle each day, but on occasion a 
single source was observed at multiple hour angles 



on the same day. The flux density scale was set 
by observations of 3C 286. 

The typical rms errors in flux density (L01) are, 
before 1989 August: 

al = 0.037 2 +(0.014S 2 ) 2 of = 0.057 2 +(0.049S* 8 ) 2 

(1) 

and after 1989 August: 

al = 0.0037 2 +(0.0155 2 ) 2 of = 0.0057 2 +(0.05S' 8 ) 2 

(2) 

Here S2 and S$ are the flux densities in the S- and 
X band, respectively. 

Attempts to estimate the magnitude of system- 
atic errors suggest that these may approach 10% 
at S-band and 20% at X-band. There is also a 
systematic ripple, up to approximately 5% peak- 
peak in a few sources, introduced by the calibra- 
tion observations of 3C286. Its origin is unknown 
(F87) , but its 1-yr timescale suggests a seasonally- 
induced instrumental effect. 

Clearly discrepant data were edited out of the 
published time series by removing those data that 
deviated by more than 10a at X band or 15a at S 
band from a running boxcar mean. Nonetheless, 
as can be seen in many of the light curves (L01, 
Figure 8), there are occasional obvious one-day de- 
viations remaining in the light curves. We do not 
regard these as valid rapid changes in flux density 
and we removed them from the time series by a 
two-step editting process. In each step each data 
point was compared with a smooth model value 
obtained from a parabolic fit to a 365-day subin- 
terval centered on that sample. If it deviated by 
more than 3.5a it was rejected; here a is the rms 
deviation between the data and the model over 
that subinterval. The rejected points have some 
influence on the a and so in a second step the pro- 
cedure was repeated, rejecting a few more points. 
With two editting passes the numbers of points 
rejected from each source in each band were typ- 
ically 10-30 in 1988-1996 but none in 1979-1987. 
A few of the rejected points were from the same 
day at S and X band, which would be due to a 
common failure in the receiver system or due to a 
calibration error common to both bands. 

We analyzed the data in two sets — 46 sources 
during 1979-1987 at 2.7 GHz and 8.1 GHz and 146 
sources during 1988-1996 at 2.25 and 8.3 GHz. As 
noted above the second set started with about 18 
months at 2.7 GHz. While the frequency change 
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caused little discontinuity for most sources with 
relatively flat spectra, there were some for which 
the frequency change caused a substantial discon- 
tinuity in flux density. In such cases we scaled 
their 2.7 GHz flux densities to 2.25 GHz using the 
mean spectral index between 2.25 and 8.3 GHz. 

After the conclusion of the GBI monitoring 
program, Fey & Chariot (1997, 2000) quantified 
the compactness of many of the sources moni- 
tored. They determined the contribution of the 
sources' structures to interferometric group de- 
lays in VLBI observations and characterized the 
sources by "structure indices" ranging from 1 to 4. 
A source with a structure index of I is relatively 
compact on scales of order I mas, having 90% or 
more of its flux density in a single component, 
while a structure index of 4 indicates a relatively 
extended source, typically having two or more 
components of approximately equal flux density. 
The structure indices were listed in Table I of LOI 
and are used in our analysis below. From models 
fitted to these and other VLBI observations in S- 
band we obtained the flux density and diameter of 
the most compact source component. However, we 
re- iterate a caution made in LOI. It is well known 
that sources may show structural changes on mil- 
liarcsecond scales on timescales of decades. Since 
the structural estimates are computed from VLBI 
observations after the end of the GBI monitoring 
program, they should be taken only as indicators 
of a source's structure during the flux density mon- 
itoring. 

3. Data Analysis 
3.1. Methodology 

We motivate our analysis by considering the 
blazar 1514+197 1 . Figure 1 shows its light curve 
from 1988.3 to 1995.3. Immediately evident is 
that the 2 GHz variations are much more rapid 
(5-50 days) than the 8 GHz variations, which en- 
compass only three main peaks in 7 yrs. These 
slow variations also appear at 2 GHz, but with 
a reduced amplitude and significant delay com- 
pared to those at 8 GHz, as expected for intrin- 
sic synchrotron bursts. In contrast, the fast vari- 
ations are relatively strong at 2 GHz and much 

Throughout the paper we refer to sources by their 1950 
coordinates 



weaker at 8 GHz, which is consistent with the ex- 
pected decrease in ISS amplitude with increasing 
frequency. In most of the sources monitored there 
are 2 GHz variations present on timescales from 
5-50 days but at somewhat lower amplitudes than 
for 1514+197. Consequently, we adopt the hypoth- 
esis that the rapid variations visible in most of 
the sources at 2 GHz are dominated by ISS and 
the slower (year-long) variations at both frequen- 
cies are intrinsic. Figure 2 shows a source at 
a lower Galactic latitude (-11°) with similar be- 
haviour but with a longer timescale at 2 GHz than 
for 1514+197, which is at latitude 56°. 

We separated the fast and slow variations by 
applying low-pass and high-pass filters to the light 
curves. For each time series we smoothed the data 
by a cosine function truncated at its first zeroes 
with a full width T sm . This smoothed series char- 
acterizes the slow variations and was subtracted 
from the raw series, leaving the difference as the 
high-pass time series. In the figures the smoothed 
data are shown by the solid lines and the high- 
pass time series is simply the deviation of the in- 
dividual points from the smooth curve. Wc ini- 
tially set T sm to 365 days at both 2 and 8 GHz, 
since typical time-scales for cm-wavelength vari- 
ations have been a year or more (e.g., Medd et 
al. 1972; Hughes, Aller & Aller 1992). For a 
few sources there were substantial slow variations 
which were nevertheless faster than the 365-day 
smoothed curve. For these sources we adopted a 
shorter smoothing time in the range 75-200 days 
as listed in Table 1. Several of these sources are 
well-known variables at cm-wavelengths and the 
separation of intrinsic variations from ISS vari- 
ations is less clear than for those with intrinsic 
times-scales of a year or more. 

In order to study the scintillation phenomenon 
in the short term variations, we want to normal- 
ize them by the mean intrinsic flux density of each 
source. Since many of the sources do vary intrinsi- 
cally, as characterized by the low-pass time series, 
we divided the high-pass series by the low-pass se- 
ries for each source. We refer to the result as the 
scintillation (ISS) data, from which we estimated 
a timescale and an rms amplitude. For each source 
we computed the auto-correlation function (acf ) at 
both frequencies and the cross-correlation function 
(ccf) between them. Though the data were sam- 
pled at nominal intervals of cither 1 or 2 (sidereal) 
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Fig. 1. Flux density at 2 and 8 GHz for 1514+197. The solid curve is the time series in Jansky smoothed 
over 365 days with a truncated cosine function. It is interpreted to be dominated by intrinsic variation in 
the synchroton radiation from the source. The short term deviations from the line are due to ISS. The trace 
below each curve gives the rms error in a single measurement (on the same scale) relative to the horizontal 
axis. 
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Fig. 2. — Flux density at 2 and 8 GHz for 0528+134 in the same format as Figure 1. 
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days, the actual sampling is not exactly regular. 
Thus we computed the acf by summing the lagged 
products into bins at lags of an integer times the 
nominal sample interval, normalizing the sum by 
the actual number of products in each bin. Fig- 
ure 3 shows examples. 

In general the acfs have a narrow spike at the 
origin, due to noise, a broader component decay- 
ing toward zero over times of 5-50 days but typ- 
ically oscillating randomly about zero at larger 
lags; superposed on this are noise variations that 
are independent from one lag to the next. It is 
the "signal" with a timescale of 5-50 days that 
we attribute to ISS. We reduced the effect of the 
noise by smoothing each acf by a Gaussian func- 
tion with a 3.3 day full width at half-maximum, 
but first replaced the zero lag point by the acf at 
the first lag (1 or 2 days). The signal variance was 
then estimated from this smoothed acf at zero lag. 
This procedure assumes that the ISS timescale is 
several days; any rapid (intra-day) variations will 
not be included in such an estimate but will con- 
tribute to the "noise" . The rms amplitude is taken 
as the square root of this signal variance giving the 
short-term modulation index m = S lms / (S), since 
we have already divided by the (varying) source 
flux density. The signal timescale r was found as 
the lag where the smoothed acf crosses half of the 
signal variance. 

This analysis, in principle, allows us to estimate 
S lms to values well below equation (2), which ap- 
plies to a single observation of flux density. How- 
ever, the limit in the equation near m = 0.015 is 
mostly due to uncertainties in the calibration of 
flux density, which is subject to low amplitude er- 
rors - some of them systematic. These may be in- 
fluenced by weather, equipment changes and other 
factors that arc not well monitored. Thus, even 
though our m estimates are corrected for the in- 
fluence of noise and other measurement errors that 
are independent from one day to the next, we have 
little confidence in identifying ISS where m is less 
than about 0.01. 

In analyzing the same data, L01 characterized 
the short-term variabilty by the logarithmic slope 
of the intensity stucture function computed be- 
tween lags of 2 and 32 days. They found values 
between and 1.2, with somewhat higher struc- 
ture function slopes at lower Galactic latitudes, 
suggesting that variations on times of 32 days or 



shorter were due to ISS. We have plotted this pa- 
rameter against our modulation index m 2 , and 
find them to be highly correlated. In general a 
value of 0.67 is expected for the structure func- 
tion slope for pure ISS in a Kolmogorov scatter- 
ing medium that extends up to the observer (see 
Coles 1988, Blandford, Narayan & Romani, 1986). 
However, the structure function estimated by L01 
was not corrected for the influence of system noise. 
Since the structure function for pure white noise 
is flat, the slope parameter will be significantly re- 
duced when the variance due to ISS is comparable 
to that due to system noise. We use our m 2 esti- 
mate to characterize the ISS, since it is corrected 
for noise. We return to the question of the struc- 
ture function slope in §5.5. 

Table 1 tabulates the results with 1950 coor- 
dinates as the source name with Galactic coordi- 
nates and effective length through the interstel- 
lar plasma (see §4.1). We omit SS433 (1909+048) 
from the Table and include the Galactic center 
source Sgr A* (1742—289), but neither are con- 
sidered in the statistical analyses that follow. We 
give mean flux densities, modulation index and 
timescale as m 2 , m 8 , t 2 and r 8 at 2 and 8 GHz. 
The formal errors in these parameters are typically 
10-20%, but the acf plots sometimes showed large 
fluctuations versus time lag. By inspecting each 
plot we assigned an "ISS confidence level" (0,1,2) 
at each frequency: 2 is for convincing evidence 
for short-term variations (errors 10-20%); 1 is for 
just credible evidence (errors 20-50%); when no 
credible short-term variations are detected (even 
though m and r were formally estimated). This 
ISS confidence level is given in brackets following 
the enries for m 2 , m&- While at 2 GHz 121 sources 
show a credible ISS "signal" and a timescale in the 
range 5-50 days, only 60 sources show credible ISS 
component at 8 GHz. For the sources observed 
in 1979-87 as well as in 1988-96, the parameters 
were estimated separately from each data set. We 
find a strong correlation (85%) between the two 
independent estimates of m 2 and somewhat lower 
correlation (42%) between the two estimates of r 2 . 
In Table 1 we list the averages from the two epochs 
for the 38 sources common to both data sets, and 
we omit the 5 sources in the first but not in the 
second data set. 

The cefs were computed in the same fashion 
but normalized by the square root of the product 
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Fig. 3. — Correlation functions for the high-pass filtered time series for 1514+197 (left) and 0528+134 (right) 
at 2 and 8 GHz. Top: Normalized auto-correlations with the 8 GHz acf plotted against negative time lags. 
The overplottcd lines are smoothed as described in the text. The signal variance is found from the zero lag 
value of the smoothed curve before normalization. The timescale estimates are indicated by vertical arrows. 
Bottom: The cross correlation of the 2 and 8 GHz time series, normalized by the square root of the product 
of the two signal variances. 



of the signal variances at the two frequencies. If 
the varying signal were fully correlated between 2 
and 8 GHz, the two acfs would have the same 
shape (though perhaps different variances) and 
the ccf would show the same symmetrical shape 
peaking near unity at zero lag. Few of the ob- 
served ccfs follow such ideal behaviours but typ- 
ically show random-looking variations over time 
lags of 10 days or more with superimposed noisy 
(day-day) variations. Several of the sources with 
large amplitude slow variations show partial leak- 
age of the intrinsic variation through the high-pass 
filter; this appears as a cross-correlation peaking 
at positive lags, which corresponds to 8 GHz vari- 
ations preceding those at 2 GHz. 

It was found also that several sources exhib- 
ited a low amplitude spike in the ccf at zero lag. 
Since these were mostly sources with an ISS con- 
fidence level of and 1, the spike is not due to 
very rapid wide-band ISS. Of course, the additive 
noise in the time series should be entirely indepen- 
dent at the two frequencies. Thus the spike is an 



unexpected correlation introduced perhaps in the 
calibration process, due to the influence of a tem- 
perature change on system gain, bad weather, or 
some similar systematic effect, which varies inde- 
pendently from day to day. To avoid this influenc- 
ing our cross-correlation estimates we smoothed 
the ccfs (ignoring the spikes) in the same fashion 
as for the acfs and recorded the smoothed ccf at 
zero lag as a measure of the cross-correlation coef- 
ficient. In the table we list this zero-lag correlation 
coefficient only for the 36 sources with ISS confi- 
dence level of two at 2 GHz and one or more at 
8 GHz. Among these several cross correlations are 
influenced by both ISS and relatively rapid intrin- 
sic variations which typically shift the correlation 
peak to positive time lags, as noted above. 

3.2. Overview of Variations 

We have completed the analysis described 
above on all of the sources and have tried to clas- 
sify their variations, by examining all of their 2 
and 8 GHz light curves (see Figure 8 of L01) and 
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the correlation plots as in Figure 3. In addition 
to listing the parameters of the short-term varia- 
tions and the ISS confidence levels, we also list in 
Table 1 a characterization of the slower variations 
as M 2 and M§, which give the rms variation on 
a 300-day timescale, normalized by the mean flux 
density. (This quantity was obtained from the 
square root of one half of the structure functions 
of the low-pass time series at a time lag of 300 
days for both 2 and 8 GHz). 

We find that there are a several types of light 
curve, that depend on the relative importance of 
ISS and intrinsic variations. For some sources 
(e.g., 0235+164, 1749+096, 2200+420) intrinsic 
variations dominate, showing a sequence of out- 
bursts on times of a few months to 2 yrs, of- 
ten delayed and reduced in amplitude at 2 GHz 
relative to 8 GHz, as expected for intrinsic syn- 
chrotron flares. For some of these with faster out- 
bursts we needed to reduce the smoothing time 
T sm (as listed in Table 1), and so the separa- 
tion of ISS from intrinsic variations is less clear. 
At 2 GHz many sources show both the long- 
term (intrinsic) variations and short-term (ISS) 
variations over 5-50 days. In a few cases (e.g., 
1635—035) there is a larger overall amplitude vari- 
ation at 2 GHz than at 8 GHz. The 2 GHz 
short-term timescale r 2 was somewhat longer for 
sources at lower Galactic latitudes (e.g. 0528+134, 
0653-033, 1741-038), but timescales of 20 days 
were also seen in a few sources at mid-latitudes 
(e.g. 0537-158, 1614+051, 1635-035). A few 
sources show the systematic annual variation (e.g., 
5% peak-peak for 0952+179), which is also seen 
at a low level (1.5% peak-peak) in the flux density 
calibrators 1328+254 and 1328+307. 

There has already been considerable analysis of 
the discrete events known as extreme scattering 
events (ESEs) in the GBI data (see F87, Fiedler 
et al., 1994, L01). The complete set of time se- 
ries in L01 show that there are occasionally other 
events or evidence of non-stationary variation at 
2 GHz. Examples are 1538+149 (see Figure 16) 
which from 1989.9 until 1991.6 showed a dramatic 
increase in variation on month-long timescales, 
and 1502+106 and 2059+034, which both showed 
several large amplitude excursions which did not 
meet the criteria for an ESE. The fact that these 
"events" were much stronger at 2 GHz than at 
8 GHz suggests an origin in changing interstellar 



propagation conditions, but we do not pursue their 
origin further. 

4. Modeling Interstellar Scintillations and 
Intrinsic Source Structure 

We now develop a model for the expected ISS of 
the sources in the GBI monitoring program. There 
are two aspects to this modeling, the ISS itself and 
the intrinsic source structure. We discuss each in 
turn, but put the mathematical details of the ISS 
in Appendix A. 

4.1. Interstellar Scintillation Model 

From pulsar observations it is known that the 
ionized ISM has a disc-like Galactic distribu- 
tion. The electron density microstructure is well- 
modeled by a power-law power spectrum with a 
Kolmogorov spectral index on many lines of sight 
(Armstrong, Rickett, & Spangler 1995). Tay- 
lor & Cordes (1993, hereafter TC93) developed 
a detailed model for the electron distribution in 
the ISM, which has been revised by Cordes and 
Lazio (2005, hereafter CL05). Both papers give an 
explicit model for the local mean electron density 
and the strength of its small scale variability as 
characterized by the strength parameter C%. This 
parameter sets the level of the density spectrum, 
P3N = C%K~ n / 3 , at three-dimensional wavenum- 
ber k. In this expression the inhomogeneities are 
described as isotropic, since along a given line of 
sight any anisotropics in the density structures 
are expected to be randomly oriented, so that the 
cumulated effect is approximately isotropic in the 
transverse plane. A measure of the cumulated 
scattering expected for a particular radio source 
is given by the scattering measure defined as the 
line of sight integral SM = / C%dl. Though the 
spectral index of the density power spectrum is 
near —11/3 for many lines of sight, Cfj and SM 
can vary by orders of magnitude over angular 
scales of a few degrees. 

There are several observable scattering effects, 
of which we are concerned with the variability in 
apparent radio flux density (scintillation). Scintil- 
lation is classified as weak or strong according to 
the scintillation index (m p t = rms/mean flux den- 
sity) for a point source, which depends on the ob- 
serving frequency as well as on the scattering mea- 
sure. If TOp t > 1, the scintillations are strong; if 
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fnpt *C 1, the scintillations are weak. The inverse- 
frequency dependence of the refractive index in a 
plasma leads to a decrease in the strength of scat- 
tering with increasing frequency. 

Interstellar scintillation was first observed in 
pulsars at relatively low radio frequencies (0.1- 
1 GHz) and was mostly in the strong regime, in 
which there are two branches of the phenomenon, 
diffractive and refractive scintillation (DISS and 
RISS as already referred to). The extremely 
small diameters of the emitting regions in pulsars 
ensures that they scintillate essentially as point 
sources showing both branches. See reviews by 
Rickctt (1990) and Narayan (1992) for discussions 
of the physics of DISS and RISS. Given the TC93 
distribution Walker (1998, 2001) calculated the 
transition frequency separating weak from strong 
scattering for extragalactic sources; he found val- 
ues ranging from about 6 GHz toward the Galactic 
poles to about 2 GHz for directions toward the in- 
ner Galactic plane. 

In contrast, extragalactic sources are observed 
typically at higher frequencies (2-15 GHz) and 
in addition their angular extent is much larger 
than that of pulsars. As a result extragalactic 
sources are often observed under weak ISS condi- 
tions; however, ISS is usually not detectable, be- 
cause their angular extent is usually sufficient to 
"quench" the scintillations, in the same way that 
the angular extent of planets quenches their op- 
tical twinkling in the Earth's atmosphere. These 
effects explain the rarity of ISS observed for extra- 
galactic radio sources. Low-frequency variability 
below 1 GHz (§1) is due to strong RISS of com- 
pact radio sources (e.g. Rickctt, 1986) while IDV 
is weak ISS (WISS). 

The spatial character of the scintillations in a 
transverse plane can be characterized by the auto- 
correlation of the intensity (fourth moment of the 
field) or from its two-dimensional Fourier trans- 
form versus transverse wavenumber (the intensity 
spectrum) . From the theory of wave propagation 
in random media (e.g., Ishimaru 1978), we can ob- 
tain a theoretical expression for the intensity sec- 
ond moment, but its evaluation is difficult under 
the most general conditions. However, the expan- 
sion for low wavenumbers leads to the solution for 
both weak ISS (WISS) and for RISS in strong scat- 
tering; the DISS solution is obtained in the high 
wavenumber limit (see Coles et al. 1987, hereafter 



CFRC). When the source has a significant angular 
extent the fine-scale (high wavenumber) fluctua- 
tions of the pattern are smeared out by the effects 
of overlapping patterns from neighboring source 
components. This explains why no DISS has been 
observed from any extragalactic sources (Denni- 
son & Condon 1981). It also means that we only 
need the low- wavenumber expansion to describe 
both WISS and RISS. In this limit the intensity 
spectrum is given by a line of sight integral of the 
product of the local density spectrum, the Fresnel 
filter that suppresses the power at low wavenum- 
bers and a low-pass cut-off due to the visibility of 
the source (CFRC, equation 11). 

Rickett et al. (1995) used this method, to in- 
terpret the ISS from IDV source 0917+624 ob- 
served at frequencies both above and below the 
weak/strong transition frequency (see also the 
analysis of IDV from quasar 0405—385, Rickett 
et al. 2002). Here we follow the same method but 
generalize the theory to an arbitrary distribution 
of scattering along the line of sight. In Appendix A 
we give details of the appropriate expression for 
the spatial correlation function of intensity as an 
integral over the distance z from the Earth includ- 
ing C 2 N {z) from the TC93 or CL05 models. 

The equations as given in Appendix A yield the 
scintillation index m; ss and spatial scale s; ss , un- 
der an assumed form for the source visibility func- 
tion as described in the next section. We convert 
the spatial scale to a temporal scale by dividing 
by an assumed velocity for the scintillation pat- 
tern. There are a few examples of IDV in which 
the relevant velocity has been shown to be consis- 
tent with the speed of the Earth relative to the lo- 
cal standard of rest (LSR) . Consequently we used 
the annual average of this LSR speed projected 
transverse to each line of sight as the model ve- 
locity ( 20-30 km s _1 ), and so obtained model 
time-scales t iss for each source. As Appendix A 
shows, the contribution of the scattering medium 
to the ISS is controlled by the effective scattering 
length L c s and the scattering angle 6 S . In devel- 
oping the necessary code for the calculations we 
made use of the subroutines described by TC93 
and CL05. 

4.2. Intrinsic Source Structure Model 

Having specified a model for the scattering 
medium, we must also specify a model for the 



10 



(U O 

xi o 
o 

s 



o S E ' > 2.00 


1 


_l_ 2.00 > S 2 > 0.50 


— □ 


□ 0.50 > S 2 


□ □ □ 


° □ 

□ □ 

- + ++ + :^\ 

+ + O o 

4- o 

+ 1 ° 

+ + + o ++ ° 


-p° ° 

+ + □ □ 
" ° +++ ° + + +++ + + 


1 






-50 50 

Galactic Latitude (deg) 













□ o _|_n 


b> 50.00 






s +i ^ + □ + 5000 > 


b> 35.00 








b 






° D □ □ 








B - » + . n 








° N + +D + + □ □ 

o ov -ff U , + □ 
































+ □ 

° X + 


□ 

□ 






N- 
















N . + 










o + 
















+ 






i 







examine three source models with peak brightness 
temperatures of 10 11 , 10 12 , and 10 13 K. For each 
of the 146 sources, we assumed 50% of their flux 
density to come from a Gaussian brightness core 
of these three temperatures. Hence we obtained 
three diameters (full width at half maximum 
FWHM) using 



'FWHM 



(l.l// G Hz)/Sjy/T b ,i 



2 mas 



(3) 



where the core flux density is Sj y , the peak bright- 
ness temperature is T bil2 in units of 10 12 K and 
#fwhm is in mas. 

Figure 4 shows the computed mi ss at 2.25 GHz 
for the 10 12 K source models under the CL05 
model for the electron distribution versus Galac- 
tic latitude and total mean source flux density (S) . 
There is only a weak trend for increased to; ss at low 
latitudes; this illustrates its dependence on prop- 
erties of both the ISM and the source. The right 
panel shows a lower envelope to the (inverse) rela- 
tion of m to (S) . It is approximated by the dashed 
line, whose slope is from CFRC (equation 13), 
which with 9 so cx y/JS) gives mi SS cx (S) 7 ^ 12 . 



Fig. 4. — Model scintillation index at 2.25 GHz 
computed under the CL05 model for the Galac- 
tic electrons, assuming that each source has 50% 
of its flux density in a Gaussian core of 10 12 K 
brightness. Model scintillation index versus Left: 
Galactic latitude and Right: mean flux density 
with dashed line mj ss oc (S)~ 7 ^ 12 . 



visibility functions for the sources in the survey. 
Rcadhcad (1994) has argued that the distribu- 
tion of brightness temperatures from relativisti- 
cally beamed jets of 5 GHz-selected radio sources 
ranges from ~ 3 x 10 10 to 3 x 10 12 K. Accord- 
ingly, we adopt a simple model in which the most 
compact component of each source has a circu- 
lar Gaussian brightness distribution in which the 
peak brightness is limited to some maximum due, 
for example, to synchrotron self absorption or in- 
verse Compton losses. Brightness temperatures up 
to about 20 times greater than the equi-partition 
value of 3 x 10 11 K are seen in VLBI images of 
many radio-loud active galactic nuclei. These are 
well described by relativistic beaming from narrow 
jets with Doppler factors up to about 20. Thus we 



5. 2 GHz Observations Compared with 
ISS Model 

In this section we confront the model with the 
observations. We begin with global comparisons 
of the observations with ISS predictions from the 
model. Next we compare the fast variations to 
the slow variations, since both the fast ISS and 
slow intrinsic variations are indicators of compact 
structure. Finally, in a limited number of cases 
where the ratio of the ISS signal to the noise is 
sufficiently strong, we perform two further analy- 
ses. First we do a power spectrum analysis and 
compare with theory. Second we epoch-average 
the timescale separately for each of the 12 months 
in the year, in a search for an annual modulation 
in the ISS timscale due to the changing Earth ve- 
locity as it orbits the Sun relative to the scattering 
medium. 

5.1. A Global Comparison for Scintillation 
Index 

Since the short-term variations are much 
stronger at 2 than at 8 GHz, we compare the 
modulation indices TO2 and timescales ti with the 
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Fig. 5. — Scintillation index at 2 GHz against mean spectral index between 2 and 8 GHz. Left: The different 
symbols flag ranges of the magnitude of the Galactic latitude. Right: Subset of sources with structure index 
1-4, as an indicator of increasing source size. 



ISS models. In making the comparison, we must 
consider the two major influences, which are the 
source structure and the distribution of the scat- 
tering plasma, for which we use the CL05 model. 
With our simplistic source model (i.e., 10 12 K 
component with 50% of the flux density), the re- 
sults are only meaningful in a statistical sense. In 
addition the CL05 model only approximates the 
Galactic distribution of the warm ionized medium 
and substantial deviations are to be expected on 
particular lines of sight. In many of the figures 
that follow we plot the results against a param- 
eter of the scattering model and flag the points 
according to some measure of the source model, 
or vice- versa. 

Figure 5 plots the observed 777,2 against mean 
spectral index (a) between 2 and 8 GHz (S oc v a ), 
which is also listed in Table 1. We see that for rel- 
atively flat spectrum sources (say a > —0.4) there 
is a much higher fraction of sources with scintilla- 
tion indices above 1%, than for steeper spectrum 
sources, as expected if steeper spectrum sources 
have larger diameters. In the left panel the points 
are flagged by the magnitude of the Galactic lati- 
tude (|6|), showing the expected trend for greater 
ISS at lower latitudes. For the steepest spectra 
(a < —0.6) there is little or no significant scintilla- 
tion, and so in the analysis that follows we exclude 
the 17 steep spectrum sources (a < —0.6), which 
include the flux density calibrators. 

The right panel shows the subset of sources with 



VLBI-observed structure index. The sources with 
structure indices of 1 or 2 have mostly flatter spec- 
tra than those with more extended structure. Also 
many of the compact structure sources show scin- 
tillation indices above 2%. While a structure in- 
dex of 1 or 2 does imply a small diameter, we note 
that structure index is defined as a measure the 
simplicity of a source's fine structure, that is par- 
tially resolved in VLBI; larger structure indices 
indicate complex structures which can confuse the 
use of VLBI observations for geodetic purposes. 
Thus it is not a simple measure of source angular 
diameter nor of the flux density in a source that is 
unresolved in a VLBI observation. We have also 
examined a similar plot flagged by the VLBI model 
source diameters mentioned in section 2, but there 
is not a clear trend of m 2 with these diameters. 

Figure 6 compares the observed scintillation in- 
dices directly with those predicted from the 10 11 K 
and 10 12 K models. In this and subsequent fig- 
ures, except where noted, we include only those 
sources with 2 GHz ISS confidence levels of 1 and 
2 and spectral index a > —0.6. The left panel 
shows that the 10 11 K model correlates weakly 
with the observations, with observed 777,2 values ly- 
ing mostly above the predictions. Note that more 
sources with flat or inverted spectra tend to lie on 
the high side of the prediction, as expected. In 
the right panel the observed m 2 values lie mostly 
below the predictions for the 10 12 K model. We 
conclude that typical brightness temperatures for 
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Fig. 6. — Observed ISS index at 2 GHz plotted against the model predictions. In the two models shown half 
of each source flux density is assumed to be a Gaussian "core" whose peak brightness is 10 11 K and 10 12 K 
on the left and right, respectively. The points are flagged by spectral index as indicated. 



the sources are in the range 10 11 K-10 12 K. 

The lack of a tight correlation is presumably 
due to departures from our simple source model. 
One cause of scatter would be a distribution in the 
fraction of flux density in the compact core. Con- 
sider the effect of a brighter core with a smaller 
diameter (following equation 3). This would in- 
crease the scintillation index, but in order to 
match an observed »S r ms and fixed parameters of 
the scattering medium, we would need a reduced 
flux density in the compact core. Quantitatively, 
an observed to 2 compatible with 50% flux density 
in a 3 x 10 11 K component could also be modelled 
as 3 x 10 12 K component with about 5% of the 
total flux density, providing that the source diam- 
eter remains large enough to partially quench the 
scintillations. For a source small enough to scin- 
tillate like a point source the lowest value for S c 
becomes equal to the observed S rms which is also 
a few percent of the total. 

The left panel of Figure 7 shows m 2 as a func- 
tion of Galactic latitude, flagged according to a. 
There is a generally similar dependence of m 2 on 
b as in the model of Figure 4. Again the influence 
of source diameter is shown by the trend for more 
flat spectrum sources to have higher m 2 values. 
We examined a similar plot for the 56 sources with 
credible values for mg but it showed only a slight 
dependence on latitude. The right panel shows 
77i 2 plotted against I a , where I a is the intensity 
of the H Q emission as observed by the Wisconsin 



Hydrogen Alpha Mapper in the one degree beam 
nearest that direction (Haffner et al., 2002). I a 
(integrated over -80 to 80 km s _1 ) is measured in 
Rayleighs and can be interpreted as proportional 
to the emission measure in the warm ionized ISM 
(1 Rayleigh corresponds to an emission measure 
of about 2.25 cm~ 6 pc at about 8000 K); thus 
it is an independent measure of the column inte- 
grated electron density squared. There is a clear 
trend of increasing m 2 with I a , and the trend is 
strongest for the most compact sources (flat or 
inverted spectra shown as squares). The dashed 
line m 2 oc VSM tx y/la gives the theoretical rela- 
tion if emission measure is proportional scattering 
measure SM. This applies if the electron density 
spectrum follows a power law with an outer scale 
which determines the density variance, which in 
turn is proportional to the square of total electron 
density (see CL05 for example). The dashed line 
is placed to approximate the apparent upper limit 
line for m 2 . 

Figure 8 shows the influence of (S) on m 2 . 
There is a deficit of points in the upper right and 
lower left of Figure 8, but we do not see the strong 
inverse relation of to 2 and (S) visible in the model 
of Figure 4. This suggests only a weak influence 
of (5*) on source diameter. Is this due to errors 
or due to inadequacies of the model? The dotted 
line is the single point error (normalized by (S)) 
from equation (2). Though in an ideal observation 
one could detect variations well below this line, as 
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Fig. 7. — Scintillation index at 2 GHz (for sources 
with ISS confidence levels 1 and 2) plotted against 
Left: Galactic latitude and Right: Ha intensity 
(I a in Rayleighs) as explained in the text. The 
points are flagged by spectral index as indicated. 
The dashed line corresponds to m oc ^/~TZ. 



noted in §3.1 we consider a lower limit for a cred- 
ible m value to be <~1%. Since the observations 
extend up to nearly 10%, most of the points are 
not much influenced by errors. The spread must 
be caused by a combination of the range of scatter- 
ing conditions, as indicated approximately by the 
Galactic latitude and the range of actual source 
brightness. The dashed line is m 2 cx (Sr 7/12 from 
figure 4, which approximates the 10 12 K model. As 
noted earlier it is likely that there is a wide range 
in the fraction of the total flux density in the com- 
pact (scintillating) components, which contributes 
to the scatter about this inverse relation of mi and 
(S). 



Fig. 8. — Scintillation index observed at 2 GHz 
against mean source flux density. The different 
symbols flag ranges of the magnitude of the Galac- 
tic latitude. The dotted line is from equation (2) 
and the dashed line is the rough lower limit for the 
10 12 K model from Figure 4. 



In Table 1 we list the flux density and FWHM 
diameter of the most compact component derived 
from model fits 2.4 GHz VLBI observations. This 
gives compact structural information on 93 of our 
sources (with diameters from 0.3 mas to 8 mas). 
We tried these values in place of the brightness- 
limited models to calculate m 2 and t 2 . Figure 
9 shows the observed mi versus this prediction, 
flagged by the FWHM 9 2 . Most of the predicted 
77i2 values are larger than observed; but for the 
diameters larger than about 2 mas the observa- 
tions approximately follow with the predictions. A 
scatter plot of the 93 VLBI diameters against flux 
density shows no inverse relation expected if the 
structures were brightness-limited. It seems that 
extracting the single most compact component in 
model fits to VLBI data is not the best way to ex- 
tract compact structural information for compari- 
son with scintillation data, since sources typically 
show structure on a range of scales. We conclude 
that the 2.4 GHz VLBI observations are of limited 
utility, except for the more resolved sources. 

5.2. A Global Comparison for scintillation 
timescale 

We now consider the 2 GHz ISS timescale t 2 . 
Figure 10 shows observed Ti as a function of b in 
the left panel. Whereas there is clear increase in 
n at low latitudes, the increase is much less than 
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Fig. 9. — Observed 777.2 against predicted from the 
ISS model using source data from S-band VLBI. 
The points, flagged according to the source diam- 
eters, only show agreement for the larger sources. 



for the models, as shown in the right panel which 
plots the observed against the predicted timescale. 
At latitudes below about 10 degrees the model 
timescales rise to more than 500 d, due to the very 
long scattering length through the inner Galaxy. 
The timescale t; ss is given by the spatial scale Si ss 
divided by the Earth's velocity with respect to the 
scattering plasma. In our model we took the ve- 
locity to be the yearly average of the Earth's trans- 
verse speed with respect to the LSR. Whereas this 
may be reasonable for scattering paths out to a few 
hundred parsecs, it will systematically underesti- 
mate the velocity at greater distances, due to the 
absence of differential Galactic rotation and the 
peculiar motions of the ionized ISM in the model. 
This presumably explains at least part of the dis- 
crepancy. 

Overall, we conclude that the general behaviour 
of the scintillation index m 2 at 2 GHz is quite well 
described by the CL05 model of the interstellar 
electrons with the flat spectrum sources having a 
substantial fraction of their flux density in regions 
of 3 x 10 11 K peak brightness. Ten times higher 
brightness is possible only from sources with less 
than about 5% of their total flux density in the 
most compact core. Our analysis shows only a 
weak inverse relation of source size to flux density. 
Thus the sources must be exhibiting a range of 
peak brightness temperatures around 3 x 10 11 K, 
as could be explained by a spread in their Dopplcr 
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Fig. 10. — Left: 2 GHz ISS timescale versus Galac- 
tic latitude, flagged for the indicated ranges a. At 
a given latitude we expect shorter timescales from 
the more compact sources which typically have 
higher values of a. Right: 2 GHz ISS timescale 
versus model r 2 for the 10 12 K source model, 
flagged for the Galactic latitude as indicated. 

factors. 

5.3. Comparison of Fast and Slow Varia- 
tions 

As already noted there are both intrinsic and 
ISS variations in many of the GBI time series. In- 
deed the presence of intrinsic (slower) variations 
was part of the selection criteria for the sample. 
Since slow intrinsic variations and fast ISS varia- 
tions are both indicators of compact source struc- 
ture, we examine the relationship between them. 

As described in §3.1 we characterize the low- 
pass filtered time series by M 2 and M 8 , which are 
the rms variations on a timescale of 300 days (nor- 
malized by the mean flux density) , as given in Ta- 
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Fig. 11. — 2 GHz comparison of short-term and 
300 day modulation indices observed at 2 GHz, 
flagged by 300 day modulation indices at 8 GHz. 

ble 1. We note that these quantities respond to 
variations on scales of 300-days or longer, includ- 
ing linear trends across the entire data span. Fig- 
ure 11 shows the 2 GHz ISS index mi against M2, 
flagged by M 8 . In this plot we included all the 
146 sources, to avoid a bias against low TO2 values 
for cases where the 2 GHz ISS was not reliably es- 
timated (ISS confidence level 0). The plot shows 
a strong correlation of the incidence of ISS and 
intrinsic variations at 2 GHz. This is consistent 
with both ISS and intrinsic variations being as- 
sociated with small diameters. The diagram also 
shows intrinsic variations to be strongly correlated 
between 2 and 8 GHz as expected. 

We also note that Hughes, Aller, & Aller (1992) 
found time-scales and slopes for the cm- wavelength 
variation from a set of sources monitored with the 
26 m antenna at the University of Michigan. 27 
of their sources are in the GBI source list. Their 
time-scales range from 0.3 to over 3 years, and 
their structure functions on these time-scales were 
clustered about 1 year, which they ascribe to in- 
trinsic variation. 

5.4. Annual Modulation in ISS time-scale 

It has been established that there are an- 
nual modulations of the characteristic timescalc 
in several IDV sources (J1819+385, 0917+624, 
1257-326 1519-273: Dennett-Thorpe and de 
Bruyn, 2003; Rickett et al., 2001; Jauncey and 
Macquart, 2001: Bignall et al., 2003; Jauncey et 



al., 2003). 

The modulations are due to the annual cycle 
of the Earth's velocity relative to the scattering 
medium. Accordingly, we have looked for similar 
annual modulations in the ISS signal at 2 GHz. 
The right panels of Figure 12 show the predicted 
annual velocity variation for blazar 1514+197. 

Our analysis consisted of computing structure 
functions out to lags of 30 days averaged for each 
calendar month over all years available. These 12 
estimates and their errors are shown by the points 
in the left panel of Figure 12 for 1514+197. (De- 
cember is repeated as month and 12). From 
plots such as this we estimated 12 timescales t a , 
by fitting smooth models to each monthly struc- 
ture funtion. Wc first fitted a form to the structure 
function of the entire time series and selected the 
best fit as defining the functional form to be used 
for each month. We tried various forms including 
an exponential and functions of the form: 

D m (r)=Ds^. (4) 

We optimized the parameters D s , t and p (in the 
range 1.4 to 2.0) and also the exponential (used in 
place of p=l) and selected the best fit. With D s 
and p then fixed, we found the best-fitting value 
of the time scale t for each of the 12 months for 
that source. When available we did independent 
analyses for the 1979-87 and 1988-96 data sets. 

Figure 13 shows some results as nine plots for 
six sources. The best example is 1514+197 in the 
lower right panel. The points are t and their er- 
rors; the lines are predictions (cx 1/Vlsr), assum- 
ing that the scattering medium is stationary in the 
LSR, where the constant of proportionality was set 
to give one half of the maximum observed t a for 
the annual average of Vlsr- For 1514+197 there is 
good agreement between the LSR-model and the 
observations. See also the right panel of Figure 
12. However, for the great majority of the other 
sources there is not a convincing trend of t a over 
the year (relative to the errors); even though, with 
many years of data, we have as many as 7-8 inde- 
pendent month-long samples for over 100 sources. 

There are several factors affecting the analysis. 
First, the reduced signal-to- noise in one twelfth 
of the data; second, when r 2 > 20 days there is 
only one independent "scintlc" in each month; 
third, for low latitude sources the long scatter- 
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Fig. 12. — Left: Monthly 2 GHz structure function estimates for 1514+197. December is repeated as month 
and 12. Right: Earth velocity relative to LSR perpendicular to the line of sight toward 1514+197 (Vlsr) 
against day of year. 



ing path makes for large differential motion (sys- 
tematic and random) through the Galaxy, which 
will tend to dominate the ~30 km s" 1 variations 
of the Earth; fourth, for several sources there are 
large changes in the flux density over one-month 
timescales which are, presumably, intrinsic and 
may correspond to changing fine structure in the 
source on similar timescales. All of these factors 
will confuse the systematic annual variation. Fur- 
ther the phase (and shape) of the annual curve 
depend on both the ecliptic coordinates of the 
source and on the LSR velocity of the scatter- 
ing plasma. Though the IDV in sources 0917+624 
and 1257—326 were well described by the LSR- 
model, the IDV of J1819+385 requires a velocity 
of about 20 km s _1 relative to the LSR. Thus our 
predicted curves are only a first-order model for 
the ISS timescale. 

Panels 1-6 of the Figure 13 show three sources 
(0133+478, 0333+321 and 0420-014) observed 
in both 1979-87 and 1988-96. They were se- 
lected as the sources observed in both periods with 
77i2 ^ 0.03 and in which there was some appear- 
ance of an annual modulation in agreement with 



the predictions. Comparing the upper and cen- 
tral panels the data points show some relation to 
the predictions in both epochs for each source. 
However, individual months have significant differ- 
ences between the epochs. The lower panels show 
the sources 1528+134, 1653-033 and 1514+197, 
which are the best examples of annual cycles seen 
in the 1988-96 data set. 

5.5. The ISM model 

In the foregoing we have shown that the rapid 
2 GHz flux density variations in the GBI data 
are generally consistent with predictions from the 
CL05 model for the electron distribution in our 
Galaxy, assuming models for the most compact 
components in the 146 extragalactic radio sources 
that do not require extreme Doppler boosting fac- 
tors and are in general consistent with the max- 
imum brightness obtained from VLBI. We have 
also compared the observations with predictions 
from the TC93 electron model. The predicted 
scintillation levels are somewhat higher from the 
TC93 model, such that the typical implied source 
brightness is 10 11 K, i.e. a factor three below that 
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RISS using NE2001 model l= 180.0 b= 45.0 



Fig. 13. — Monthly 2 GHz timescale estimates for 
six sources: top panels from 1979-87; center and 
bottom panels from 1988-96. Sources 0133+478, 
0333+321 and 0420-014 were observed in both 
periods. The points show observations; the line 
shows the LSR prediction. 

implied using the CL05 model. This difference 
highlights the importance of refining our knowl- 
edge of the electron distribution in the estimation 
of source structure from ISS. 

As mentioned in the introduction variation on 
timescales as short as one hour are seen from 
a few extragalactic sources at cm wavelengths. 
Such rapid variations can only be explained by lo- 
cal scattering "clouds" located 10-30 pc from the 
Earth. This is in marked contrast to our results 
reported here in which the CL05 model can suc- 
cessfully explain the 2 GHz variations from 121 
comapct radio sources. The CL05 model involves 
scattering at typical distances of 500-1000 pc and 
has a deficit of electrons in a "local bubble" ; thus 
it specifically excludes scattering in the local ISM. 




-0.1 1 10 

time lag (days) 

Fig. 14. — Theoretical ISS structure func- 
tions for sources with 0.5 Jansky in a Gaussian 
core of brightness 10 n K (solid), 10 12 K (dashed) 
and 10 13 K (dash-dot). Calculations are for 
180°Galactic longitude and 45°latitude with the 
CL05 electron model. The lower steeper curves in- 
clude a local bubble with much reduced electronso 
that the electron distribution extends to the solar 
vicinity. 

The importance of this local bubble on the pre- 
dicted ISS is illustrated in Figure 14. Here the 
lower curves have reduced scattering plasma in 
the local bubble (i.e. within about 100 pc of the 
Sun). Related plots were made by Blandford, 
Narayan and Romani (1986) and by Coles (1988). 
The figure shows that observations of the loga- 
rithmic slope of the structure can be used to test 
for scintillation caused in the local ISM; slopes 
as steep as two are expected when the local bub- 
ble eliminates any local scattering, and slopes of 
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0.67 are expected if the scattering extends to near 
the Sun. In L01 the slopes estimated from the 
2 GHz GBI data were found to lie in the range 
to 1.2, however, these estimates were not cor- 
rected for noise and so cannot be used for this 
purpose. We do not think that there is a sufficient 
range of timcscalcs in the data to reliably deter- 
mine these slopes. The plots of structure function 
Si(T) are useful in predicting the apparent scin- 
tillation index m in flux density over an observ- 
ing span T, as m ~ ^/0.5S7(T). The differences 
in Figure 14 largely explain the higher predicted 
scintillation indices from the TC93 model, which 
does not model the local bubble. 

An important result bearing on the existence 
of local clouds comes from the MASIV survey re- 
ported by Lovell et al. (2003). They observed 
710 flat spectrum sources at 4.9 GHz for 72 hours 
and found 85 with significant variation over that 
time. However, they only detected one source with 
variations as fast as 4 hours, which was the well- 
known IDV quasar J1819+385. The sources sur- 
veyed covered most of the Northern sky and so 
suggest that the covering fraction of local scatter- 
ing clouds is small (1/85). The GBI observations 
were not sensitive to ISS on times shorter than 
one day, but the rarity of time-scales in the range 
2-4 d also suggests that the local ISM is indeed 
depleted of scattering plasma and so supports the 
local bubble concept. 

Another aspect of the three very rapid variables 
is that detailed modelling of their ISS shows that 
their local clouds cause anisotropic scattering, i.e. 
elongated plasma density structures presumably 
oriented parallel to the local magnetic field. In 
contrast we assume the scattering to be circularly 
symmetric in our modelling of the GBI data. Es- 
sentially this assumes that the length scale for ran- 
domization of the projected magnetic field is sub- 
stantilly smaller than the 500-1000 pc distances 
through the ISM. However, there is insufficient in- 
formation in the GBI data to test this assump- 
tion. Anisotropy would cause the structure func- 
tions (as in Figure 14) to have different widths and 
shapes in orthogonal directions. These considera- 
tions should be included in analysis of any rapidly 
sampled flux densities in future observations (such 
as that proposed at the PisGah Observatory). 



6. Comments on Individual Sources 

We examined plots as in Figures 1 and 3 for 
all of the 44 sources observed in 1979-87 and the 
146 sources observed in 1988-96. Whereas the pa- 
rameters derived from the correlation analysis are 
listed in Table 1 and whereas most conform to 
our hypothesis of rapid ISS at 2 GHz and slower 
intrinsic variation at 8 GHz, we now discuss some 
sources which had unusual time series or were oth- 
erwise notable. The interested reader is encour- 
aged to peruse the 146 time series plots in L01; 
where there is also a list of identifications and al- 
ternate names for the sources. 

Among the sources discussed individually in 
this section 7 out of 12 are BL Lac objects, com- 
pared with 20 BL Lacs among the 146 sources 
monitored. This may suggest that unusual ISS 
or unusual intrinsic variations are more common 
in the BL Lacs than in the quasars, but such be- 
haviour was part of the selection criterion for the 
observations. 

0016+ 731 

This quasar shows obvious variations (m 2 = 
5%) on a 1 month timescale at 2 GHz with no slow 
variations at 8 GHz. The 2 GHz variations appear 
to be relatively slow ISS; at a latitude of 11° it is 
near the upper envelope of the Ti points in Fig- 
ure 10. It has a relatively high effective scattering 
distance of 2.55 kpc and so is probably exhibiting 
refractive ISS. At 8 GHz there may also be ISS 
variations on a timescale of 11 d. 

01 33+476 

The quasar 0133+476 is an excellent example 
of rapid ISS at 2 GHz and year-long intrinsic vari- 
ations at 8 GHz. It was sampled daily in 1979-86 
and every other day in 1988-96. In both 2 GHz 
data sets there are stochastic variations by 10% on 
timescales of 10 days; these have much lower am- 
plitudes at 8 GHz but are revealed by cross cor- 
relation. The 10-day timescale of its 2 GHz ISS 
is short considering its Galactic latitude of -14°. 
The 8 GHz data show 25% variations (presumably 
intrinsic) on timescales of a year, having a much 
lower amplitude at 2 GHz. 
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0235+164 

The BL Lac object 0235+164 is a much-studied 
variable at many radio wavelengths. Both the 

2 and 8 GHz time series are dominated by deep 
(M2 ~ .4; Mg ~ 0.6) relatively slow variations, 
which we interpret as intrinsic. The major out- 
bursts over year-long timescales are often accom- 
panied by lower level events on times of a few 
months. The 2 GHz variations are delayed by 1- 

3 months, with some evidence that the delay in- 
creased between 1982 and 1994. At 2 GHz there 
are also low level rapid variations particularly near 
the peaks in the outbursts (e.g. in 1992-93). We 
presume that these latter were ISS due to the 
emergence of new components (like those in VLBI 
images) compact enough to scintillate, which then 
expanded and quenched the ISS. 

In Figure 15 we plot the flux densities overplot- 
ted with a 50-day high pass running mean. The 
middle curve is a running estimate of Ui ss , which 
is the net rms deviation in flux density from a 
quadratic curve fitted to the 50-day interval cen- 
tered on each point. At 2 GHz <T[ SS is well above 
the lower curve which shows the rms due to noise 
and estimation error, which has not been sub- 
tracted from <7i ss . Often the ISS peaks earlier in a 
burst than does the total flux density, as expected 
for the ISS of an emerging compact emission fea- 
ture. A clear example is the major burst during 
1992-1993, where <Ji ss increases abruptly in 1992.7 
just as the 2 GHz flux density starts to rise. A 
simple model for the ISS of such a single transient 
component of diameter 0(t) and flux density S&(i) 
is 

a iss ~ S b (t)m 9 /^6{ty + 0*, (5) 

where o is the critical diameter for suppression 
of the ISS and m is the scintillation index of a 
point source (~ 1.0 at 2 GHz). The physics of 
the expanding feature will determine the relation 
of 6{t) to Sb(t); while it is optically thick one may 
expect increasing Sfe(i) cx 9(t) 2 and then a decay 
in Sft(t) as it becomes optically thin. 

The 1992.7 event suggests a very rapid rise in 
Sb (t) while its diameter <0 O . We shortened the 
time for the running estimate of <7i ss to 30-d and 
made a scatter plot of <7j ss against S(t) during the 
rising part of the burst. This revealed that ai ss 
rose from 0.06 to 0.12 Jy while S(t) increased from 
1.7 to 1.77 Jy. In other words cr; ss increased by 



0.06 Jy while the flux density increased by 0.07 Jy 
(both figures have errors of about 25%). Thus 0.07 
Jy flux density emerged with a scintillation index 
of 0.06/0.07 <~ 0.9 implying a diameter <0 O . which 
we estimate as about 10 //as using the L e ff tabu- 
lated in Table 1. An emission component of 0.07 
Jy with a diameter of 10 //as at 2.25 Ghz implies a 
remarkably high brightness temperature of about 
1.7 x 10 14 K during the few days that diss showed 
the rapid rise. Subsequently <t; ss remained con- 
stant while the total flux density nearly doubled 
from 1.77 to 3.4 Jy. This suggests a relationship 
between the effective diameter of the new compo- 
nent and its contribution Sb(t) to the total flux 
density, but further theoretical work is needed to 
make a quantitative model of the expansion. 

It is notable that Kraus et al. (1999) were moni- 
toring the source at the VLA and observed a peak 
in the flux density of 0235+164 at 1.5, 4.9 and 
8.4 GHz on modified Julian days 8904-8905, coin- 
cident with the abrupt increase in <Ji ss at 2 GHz. 
They considered both intrinsic and extrinsic expla- 
nations, and our analysis suggests the emergence 
of a very bright component at that time, which 
contributed to the rapid change in the light curve 
at 1.5 to 8.4 GHz. However, the subsequent rapid 
variations of the flux density are clearly strongly 
influenced by ISS. This event highlights the im- 
portance of both phenomena in any detailed in- 
terpretation of such outbursts. 

Viewed over the full span 1979-1996, the 2 GHz 
flux density for this source shows 30-50% varia- 
tions on times of several months combined with 
only a low average level of ISS (~ 2%) variation 
on times of 3 days. This is a lower level of ISS 
than seen for several sources with very little intrin- 
sic variaion (e.g. 0133+476 and 1514+197), and it 
suggests that even though the source is highly vari- 
able the average flux density coming from sub-mas 
structure is low relative to that in sources which 
show a higher average value for mi- 

0355+508 

This low latitude (b = -1.7°) quasar (NRAO 
150) was one of two objects from the 1979-87 
observations analyzed for ISS by Dennison et al. 
(1987). At 2.7 GHz they found a slow decay over 
7 yr and 5% variations on times of 1 yr, which 
they showed is consistent with strong refractive 
ISS as expected at such a low latitude, quenched 
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Fig. 15. — Flux density vs time for 0235+164. The middle curve is a 50-day running estimate of o-; ss 
described in the text and the lower curve is the rms error in flux density; both curves scaled times 3 as 
indicated by scale on the right. 



by a 10:1 factor due to the angular diameter of 
the source which they estimated as 7 mas, which 
is comparable to the 5 mas listed in Table 1. Even 
at 8.1 GHz the model predicts refractive ISS for 
this source. 

They made no precise definition of timescalc 
but 1 yr appears to be the "quasi-period" of 
the variations about a slow decay, suggesting a 
timescale in our definition of about 6 months. 
When analyzed through our 365-day high-pass fil- 
ter we find a time scale of 1-2 months, which illus- 
trates that the use of the high-pass filter to sep- 
arate ISS and intrinsic variations is only possible 
when the two timescales are widely different. We 
detect 3% variations at 8 GHz on a timescale of 
13 days, which is comaparable to the 25 days pre- 
dicted in their model. 

0528+134 

As shown in figure 2, this quasar was quiescent 
at 8 GHz over 1988-1992. It then started a slow 
increase about doubling its flux density by 1994. 
Then over 1995.2-1996.0 it rose to 8 Jy, doubling 
its flux density again in less than 1 yr. The ISS 
at 2 GHz is 5% on a timescale of 12 d persisting 
throughout the 8 years of monitoring. We looked 
for the associated increase 2 GHz ISS amplitude 



as the 8 GHz burst began, as seen in 0235+164. 
However, though the flux density at 2 GHz did 
start to increase, the ISS amplitude did not, al- 
though it became somewhat faster. This shows 
that an increase in ISS is not always seen at the 
onset of an intrinsic burst. 

0716+714 

0716+714 is a BL Lac object studied exten- 
sively by the Bonn group and exhibits IDV at 
cm wavelengths (e.g. see Quirrenbach et al., 2000). 
Their structure function analysis shows significant 
variations on times 2 days at the 1-2% level, at 
1.4 and 4.8 GHz. In the GBI program it was 
mostly sampled at 2 day intervals and shows sig- 
nificant 2.2 GHz variation at about 4.8% rms on 
a timescale of 5 days. 

Wagner et al. (1996) also report a relationship 
between its optical and radio variations, which 
they interpret as evidence for intrinsic IDV at cm 
wavelengths. If as they suggest the 4.8 GHz IDV 
were intrinsic with a Doppler factor <~ 100, the ob- 
servable angular diameter would be ~ 3/zarcsec. 
Scaling this to 2.2 GHz, we obtain ~ 6/iarcsec, 
which should scintillate as a point source, showing 
both refractive and diffractive ISS. However, our 
2.2 GHz observations suggest refractive ISS only 
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with possible smoothing by a source diameter of 
~ 0.5 mas. 

0851+202 

Another BL Lac object 0851+202 (OJ287) is 
similar to 0235+164 in showing vigorous intrinsic 
variability on timescales of a few months, which is 
highly correlated between 2 and 8 GHz and only 
2.8% rms in ISS at 2 GHz. There are peaks in 
the short-term rms in 2 GHz flux density at the 
beginning of the intrinsic bursts, but the evidence 
is not as strong as for 0235+164. 

0954+ 658 

This BL Lac object shows 10-20% variations 
over a year at 8 and 2 GHz, with the time delay 
from high to low frequencies characteristic for in- 
trinsic bursts. It is well known since it produced 
the best and most definitive ESE in 1981 (Fiedler 
et al., 1987a). Our analysis also shows that it 
clearly exhibits typical 2 GHz ISS on a timescale 
of 4 days at 3.4% scintillation index. 

1413+135 

This BL Lac is remarkable for its strong intrin- 
sic variations (over 25%) at 8 GHz and remarkably 
low level ISS (~ 0.6%) at 2 GHz, where it also has 
only a low level of 300-day variation. Thus it ap- 
pears as an outlier near the center top of the right 
panel in Figure 11. This source has an inverted 
spectrum (a = 0.33), making a possible explana- 
tion that the compact core is heavily absorbed at 
2 GHz and contributes little to the 2 GHz flux 
density. 

1514+197 

This BL Lac object shows the highest level of 
ISS at 2 GHz (7.5% rms on a tmescale of 7 d) 
of all of the extra-galactic sources in the sample. 
At a latitude of 56° it is expected to be in weak 
ISS or near the transition frequency. Thus its high 
level of ISS suggests a smaller diameter than other 
sources. 2.3 GHz VLBI observations reported by 
Fey and Chariot (2000) list its structure index as 
1 — the most compact but its VLBI-modelled 
angular diameter in Table 1 is 1.8 mas which is 
not unusually small. The relatively high level of 
ISS in this source revealed an annual modulation 
in its ISS timescale agreeing in phase with that 



predicted for a scattering plasma at rest in the 
LSR. We have computed the power spectrum of 
its 2 GHz ISS and compared it with predictions 
for a Gaussian source model, assuming a screen or 
the CL05 and TC93 models of the ionized ISM. 
These give agreement in the width at half power 
of the ISS spectrum for a source diameter ~ 0.25 
mas, which corresponds to a brightness tempera- 
ture of <~ 10 12 K. As shown in Figure 3 and listed 
in Table 1, its ISS is also detected at 8 GHz with 
an amplitude of 2% on a timescale 11 d and is 
52% correlated with the 2 GHz ISS. This is the 
best example of 8GHz ISS in our data. 

1 538+ 149 

The 2 GHz time series in Figure 16 has an 
episode of strong variability starting abruptly in 
1989.9 and ending abruptly in 1991.6. Before and 
after this episode there are typical ISS variations 
at a few percent over 10 days; during the 20-month 
episode there were about ten excursions of ±25% 
on time scales of about one month. At 8 GHz 
there are 10% variations on times of 6 months and 
during the 2 GHz episode some signs of similar 
variations at a lower amplitude. The higher am- 
plitudes at 2 GHz strongly suggest ISS as the cause 
of the 2 GHz episode. Is this due to a change in 
the source or the medium? While we cannot give 
a definitive answer, the stochastic appearance and 
slower than normal timescale (at a latitude of 49°) 
suggests that an unusual region of enhanced tur- 
bulence intervened during 1989.9-1991.6, causing 
strong refractive variations. The 20-month dura- 
tion at 30 km s _1 implies a size of the region ~ 10 
AU. 

It is possible that this episode is similar to an 
ESE, but that the ISM structure was larger in- 
cluding multiple regions with sufficient refracting 
power to focus or defocus the 2 GHz radiation. 
Hence we looked at all sources for non-stationary 
episodes and found unusual 2 GHz variations 
in 0355+508, 0537-158, 0836+710, 1100+772, 
1502+106, 1741-038. 

1741-038 

Hjellming and Narayan (1986) analyzed one 
year of observations at 1.5, 5, 15 and 22 GHz 
of this quasar, which lies at a low latitude (13°) 
toward the inner Galaxy. They identified 20-d 
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Fig. 16. — Flux density at 2 and 8 GHz for 1538+149 in the same format as Figure 1. Note the episode of 
enhanced variation at 2 GHz from 1989.9-1991.6, as discussed in §6. 



variations at 1.5 GHz as refractive ISS (possibly 
caused at a few hundred parsecs in the North Po- 
lar Spur), while the higher frequencies were dom- 
inated by intrinsic variation. Our 2 GHz results 
also give an 18-d timescale at an rms amplitude 
of 10%. However, the 2 GHz auto-correlations are 
strongly influenced by the inclusion of data from a 
well studied ESE (1992.2-1992.6, see Lazio et al., 
2000). During the event the 2 GHz flux density de- 
creased to a minimum 50% of normal, remaining 
lower than 75% for about 2 months. 

Examination of the entire 2 GHz time series re- 
veals ±10% variations on similar timcscales from 
1983.5-1994.5. There is also a period 1988.0- 
1989.6 during which the timescale appeared to be 
much shorter (~ 12 d). The VLBI observations 
following the minimum in the light curve at 1992.4 
showed that the source had a diameter of 0.9 mas 
during the minimum flux density and 0.6 mas af- 
terwards. Lazio et al. (2000) concluded that the 
increase in diameter during the event was due to 
an extra stochastic broadening. 

For refractive ISS we expect a timescale ~ 
LQ/ 1.7V for a scattering at distance L of a source 
observed to have a diameter 9 with transverse ve- 
locity V. For 0.9 (0.6) mas at 1 kpc and 30 km s" 1 
we obtain 32 (21) d. These are within 50% of 



our estimated ISS timsecale and suggests that the 
conditions for the normal ISS and the ESE arc 
fairly similar. However, as noted for several of 
the sources discussed in this section, there is good 
evidence for non-stationary behavior in the ISS 
phenomenon, which is presumably based in dis- 
tinct physical structures crossing the line of sight 
on occasion. 

2005-403 

This quasar lies at 4 degrees latitude in the 
Cygnus region and is substantially broadened by 
interstellar scattering. As shown by Desai & Fey 
(2001) its scattering disc has a diameter of 18 
mas at 2.3 GHz with an axial ratio of (0.76:1), 
and its diameter scales as (wavlength) 2 . The as- 
sociated refractive ISS should have a time scale 
of about 1.1 years and cannot be distinguished 
from the slow intrinsic variation and (spurious) 
annual variation visible for this source. However, 
at 8 GHz there is a 4% variation on a timescale 
of ~ 12 d, which we suggest is ISS. Extrapo- 
lating their diameter measurements to 8.3 GHz 
gives a scattered diameter of maa 1.3 mas. If 
the intrinsic source size was no larger than this 
we would expect the ISS timescale to be about 
18 days xLk P c#mas/(Vkms/50), where the scatter- 
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ing path is Lkpc and Vkms is the transverse velocity 
of the medium relative to the Earth. The nom- 
inal scattering distance is 5.5 kpc, which would 
only agree with the 12 d timescale if Vkms ~ 500 
km s _1 . This is unreasonably high given a 220 
km s _1 Galactic rotation velocity, and suggests 
that the effective distance to the scattering plasma 
is much less than 5.5 kpc. At 1 kpc the timescale 
implies a transverse velocity of 90 km s -1 — still 
higher than the 30-50 km s _1 commonly assumed. 
However, these figures are necessarily crude, since 
they do not include effects of anisotropic scatter- 
ing or the distribution of scattering along the path. 
We also note that the time-series of 8 GHz flux 
density has several intervals with quasi-oscillatory 
variations which have timescales noticeably longer 
than the 12 d derived from the correlation analy- 
sis. 

7. Discussion and Conclusion 

We have used a simple filter on the time se- 
ries of flux density to separate slow intrinsic vari- 
ations from fast ISS on all of the light curves from 
the GBI monitoring program. The results confirm 
that both processes are widespread and must be 
considered together in a full description of com- 
pact radio sources. Both processes imply small 
diameters for the sources. 

Our model for the interstellar scattering is that 
it is caused by "turbulence" in the ionized ISM 
as described by the CL05 model. We also com- 
puted the ISS based on the TC93 model which 
predicted generally slightly higher scinillation in- 
dices and about 15% shorter time scales. This was 
particularly evident in the form of the predicted 
sptructure function at small time lags. 

• We find that a simple filtering can separate 
the fast ISS from the slow intrinsic variations 
in flux density at 2 GHz. A demarcation 
timescale of about 365 d is satisfactory for 
most of the sources. 

• We identify the 5-50 d variation in the 
2 GHz flux density of about 121 sources 
as ISS. These variations are the expected 
ISS phenomenon intermediate between LFV 
and IDV. 

• The modulation indices and timescales show 
clear dependence on Galactic latitude. This 



confirms that ISS is the origin of such fluctu- 
ations. In addition there is a strong correla- 
tion of the 2 GHz scintillation index with the 
emission measure estimated from WHAM 
H a emission. 

• Over the 121 sources the 2 GHz ISS mod- 
ulation index is correlated with the mean 
spectral index and also with the detection 
of compact structure on VLBI scales. 

• A model for ISS in the CL05 (or TC93) 
medium is reasonably successful if the peak 
brightness of the compact core of each source 
is in the range 10 n K to 10 12 K and 10-50% of 
the source's flux density is in such a compact 
component. This is in general agreement 
with the distribution of brightness tempera- 
tures proposed by Readhcad (1994). 

• There is a weak inverse correlation of scintil- 
lation index with mean flux density, as ex- 
pected if the effective diameter of a source 
is limited by a maximum in brightness tem- 
parature. 

• There is a significant population of sources 
that show substantial 2 GHz ISS, but do not 
vary intrinsically (above about 4% on times 
of a year or more). This implies that some 
very compact sources do not vary intrinsi- 
cally. Conversely, there are a few 8 GHz 
sources with strong intrinsic variations with 
surprisingly low levels of ISS at 2 GHz. 

• An annual modulation in the ISS time-scale 
is found for a few sources, confirming both 
the ISS model and that in these directions 
the scattering plasma moves approximately 
with the local standard of rest. However, 
most sources do not show such a simple an- 
nual cycle, implying that we do not have an 
accurate understanding of the velocity of the 
scattering plasma. 

• An intrinsic outburst from the BL Lac ob- 
ject 0235+164 was detected in 1992.7 by the 
sudden onset of ISS at 2 GHz while its to- 
tal flux density started a steep rise. We esti- 
mated the brief emergence of a component as 
bright as 1.7 x 10 14 K which decreased over a 
few days as it expanded, eventually doubling 
the total flux density of the source. 
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• At Galactic latitudes less than about 15 de- 
grees from the plane the ISS model needs to 
include differential Galactic rotation, which 
will increase the effective velocity of the ISM 
relative to the model and so decrease the pre- 
dicted time-scale, This may include a wide 
enough velocity range along the line of sight 
to invalidate the "frozen pattern" assump- 
tion. 
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A. Intensity correlation for WISS and RISS 

In the ISS model of Rickett et al. (1995) a Gaussian function was assumed to describe the dependence 
of on distance from the Earth; here we use the same basic method, but instead integrate the CL05 or 
TC93 model in the direction toward each source. 

For a source of unit flux density the intensity covariance function (at offset s) is given by CFRC equation 
(11), which we rearrange in terms of distance z measured from the observer: 



r-OC />0O 

Csi{s) = 8nr 2 e X 2 dz d 2 qC 2 N (z)q 

JQ JO 



a-2 x 



\V(qz/k)\ 2 sin 2 [^] exp[-D r (q, z) + iq ■ s] (Al) 

Here k = 2n/X for wavelength A, r e is the classical electron radius, a = 5/3 for the assumed Kolmogorov 
spectrum of density as C 2 j(z)q~ a ~ 2 at transverse wavenumber q. The refractive cut-off is given by: 

/>oo 

D r (q,z)= / dt D'(qmm(t,z)/k,t) = (qL eS (z)9 s ) a (A2) 



and D'{s,z) is given by equation (4) of CFRC. The second part of equation A2 involves an effective path 
length through the scattering medium: 

/>oo 

L cS (z) a = / dt [mm(t, z)] a C 2 N (t) /SM 
Jo 

r-OO 

SM= / dtC 2 N {t) (A3) 



It follows that near the observer the scattering distance is effectively L c g(z) ~ z; but at large distances from 
the observer L e s saturates at about the 1/e scale of the C% distribution. The scattering angle 9 S corresponds 
to the effective angular broadening at the observer: 



= k~ 



2 2,2^ r(l-a/2) 



8TT 2 riX 2 SM 



T(l + a/2)a2° 



-| -l/a 

(A4) 



We proceed as in CFRC to approximate the exponent a in A2 by a value two in place of 5/3. As this 
only governs the precise form of the refractive cut-off it changes the result very little. We further assume 
that the source visibility is a circular Gaussian function: 

V{u) = cxp[-0.5(ufc(9 so ) 2 ] (A5) 

in which the source diameter is 2.35 9 so . These two assumptions make a Gaussian cut-off in the wavenumber 
integral, which allows the integral to be done as in the Appendix of CFRC in terms of a confluent hyper- 
geometric function. We then numerically integrate equation Al with the TC93 weighting for over z to 
obtain Csi(s). 

The modulation index is then given by m — \J Cgi(0) and the spatial scale Si SS is found by solving for 
where the correlation falls to 0.5 of its maximum: Csi(s iss ) = 0.5Csi(0). 
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0.76 


1.78 


2352+495 
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